Hybrid metal-dielectric ring resonators for optical magnetic metamaterials down to 

ultraviolet range 
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In this paper, we derive a model from Maxwell equations for the magnetic resonance of ring 
resonators. Using this model we revisit the scaling of split-ring resonators. Inspired by our model, 
we propose a hybrid metal-dielectric ring resonator composed of high index dielectric material (e.g. 
SiC, Ti02, ZnS) (for the major portion) and metal (e.g. Ag). Such a new magnetic metamaterial 
is able to overcome the saturation problem of split-ring resonators and therefore is able to operate 
at short wavelength down to ultraviolet range. 
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INTRODUCTION 



Magnetic response of nature material is usually very 
weak. With the help of metamaterials, people managed 
to realize strong magnetic response from microwave fre- 
quencies to optical frequencieai"— . Optical magnetism 
can be realized through metallic structures supporting 
plasmonic modes able to interact with magnetic field, 
such as split-ring resonators (SRRs)'^~^'^ i^^'^^ and their 
derivative structures^i^iii. However, the operation of 
metallic structures at optical frequencies is limited by 
the kinetic energy of the electrons in the metal, lead- 
ing to saturation of the magnetic response when we push 
the operating frequency deeper into the optical frequency 
by size scalingi^. Although there is a sustained effort 
to push magnetic response to shorter and shorter wave- 
length, there is still lack of elegant designs for realization 
of magnetisms at short wavelengths, such as blue, violet 
and ultraviolet region. 

In this paper, we first revisit the scaling of the split- 
ring resonator and analyze the origin of the saturation of 
the magnetic response. Unlike Ref. [l^, where the kinetic 
inductance Le is derived from the kinetic energy of elec- 
trons, we derive the effective inductance from Maxwell 
equations and the dispersive permittivity of metal. From 
the revisit, we find out an effective inductance Lc which 
is the same as Le for frequencies far below the plasma 
frequency of the metal but significantly larger than 
for visible wavelength rendering scaling to short wave- 
length even more difficult. Then we propose a hybrid 
metal-dielectric ring resonator composed of high index 
dielectric material (e.g. SiC, Ti02, ZnS) (for the major 
portion) and metal (e.g. Ag). Such a new magnetic meta- 
material is able to operate at short wavelength down to 
ultraviolet range and the operating principle is inspiring 
for more designs of short wavelength metamaterials. 



II. 



MODEL DESCRIPTION 



The first artificial material exhibiting strong magnetic 
response is a type of ring resonator called split-ring 
resonator—, which sees lots of variations, with the key 
mechanism unchanged. The unit cell (a x a) of the ring 
resonator metamaterial explored in this paper is shown in 
Fig. [H It is a two dimensional ring resonator with radius 
r = 0.3a and ring thickness b = 0.2a. The ring is com- 
posed of two types of materials with permittivity ei and 
£2 , occupying 20i and 262 portion by angle, respectively. 
We assume that |e2| is large enough to ensure that the 
electric field is strongly confined inside the ring so that 
the electric field inside the ring has no radial component. 
We also assume a uniform distribution in the radial di- 
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FIG. 1. The cross-sectional schematic of the two dimensional 
ring resonator in a unit cell ax a. The radius r {r — d—^b) and 
thickness b of the ring is r = 0.3a and b = 0.2a, respectively. 
The permittivity of the 29i /202 angle part of the ring is ei /e2 . 
The background is air. 
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rection of the electric field inside the ring, so that the 
displacement current in the ring can be expressed as 



I = b——- = —iLobeiEi = — ia;6e2i?2- 
ot 



According to Maxwell equations, we have 
2reiEi + 2re2E2 = iw($'"* 



(1) 



(2) 



where = {l — F)LgI is the magnetic flux produced by 
the induced displacement current and = ^qHqFS 



{S 



is the area of the unit cell) is the external 



driving magnetic flux by incident magnetic field Hq. De- 
polarizing field contributes to the {1 — F) factor of 
where F — -kt^ ja} is the filling factor. It accounts for 
the magnetic coupling between unit cells. Then we de- 
rive from Eqs. [T]and[2]the following equation 
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+ (l - F)i-iiuLg) 



I = iu}^" 



(3) 



where C — + ^ is the effective capacitance, 
which is the serial capacitance of Ci{— (bei) / {2r9i)) and 
C2(= {be2)/{2r92)); Lg is the effective magnetic induc- 
tance of the ring, which can be approximately evaluated 
as HqFS. As long as the real part of C is positive, the ring 
will exhibit a magnetic resonance when uj = l/y^LgC. A 
magnetic resonance is possible for any of the following 
three cases: (1) ei > 0, £2 > 0; (2) ei > 0, £2 < 0; 
(3) ei < 0, e2 > 0. The dielectric magnetic metama- 
terials belong to the first case, which can be applied in 
microwave or terahertz ranges where large permittivity 
of the order several tens of eg can be easily achieved, 
but is a challenging issue in the infrared and visible fre- 
quency range due to low refractive index of materials in 
this range^i. Only recently, optical magnetism in the 
midinfrared was realized with tellurium dielectric cubic 
resonators^^. Optical magnetism in visible can also be 
realized using single crystalline silicon^^!^, but the fabri- 
cation method is quite limited in order to guarantee that 
the constituent silicon is single crystalline. The conven- 
tional metallic SRR at optical wavelength is of the second 
case, which will be studied in detail in Sec. IIIII to revisit 
the problem of saturation. The third case is a brand new 
type, which is very intriguing and significant owing to 
its ability to overcome the saturation problem of conven- 
tional SRRs for magnetic response at short wavelength 
and will be studied in detail in Sec. IIVI 



III. SCALING OF CONVENTIONAL 
SPLIT-RING RESONATOR 

To analyze the scaling of the conventional SRR at opti- 
cal wavelength, we use Ag as metal material £2 and air as 

the dielectric material ei . The permittivity of Ag can be 
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described by Drude model £2 — eo(eoo — ; 



) , which is 

dispersive and negative below plasma frequency ujp. We 




300 500 700 
Frequency (THz) 



900 



900 



„ 700 

N 

X 
H 

>-500 

d 

a> 

^300 



100 



(b) 




— g(m) 




- - g(Ki) (Lc=Le) 



100 300 500 700 900 
Frequency (THz) 



FIG. 2. (a) The comparison between the effective inductance 
Lc (solid line), the well-known kinetic inductance Lc (dash 
line) and the geometrical inductance Lg (dot line), (b) Deter- 
mination of the magnetic resonant frequency (denoted by the 
circle) with dispersive effective inductance Lc- The structure 
parameters are: a — lOOnm, r — 30nm, 6 = 20nm, 9i = 20°. 



determined the parameters £00 (= 4), Wp(= 1.4 x lO^^rad) 
and 7(= 9.69 x 10^"^rad) by fitting to the experimental 
data^^ in the frequency range 400 — 900 THz. Note that 
due to the interband transitions, the fitted Drude model 
is no longer accurate enough above 700 THz (i.e. be- 
low 430 nm wavelength), but we will still use the Drude 
model above 700 THz for principle demonstration. Since 
£2 is dispersive, C2 is also dispersive. Because of the 
negative real part and the imaginary part of £2, C2 is 
a negative effective capacitance in parallel connection 
with a resistance, which is equivalent to a positive in- 
ductance Lc in series connection with a resistance R, i.e. 
(— iwC2)^^ = —iutLc + R- Then the effective magnetic 
permeability can be expressed as a resonant form 



F'uj-' 



Moff 



(4) 



(Lg+Lc)Ci 



where F' 



F and R 



Lc is dispersive as 



shown in Fig. 2(a) (solid line). The resonant frequency 
of the effective permeability can be determined by the 
intersection point of the line g{ijj) = 1/ \J{Lg + Lc)Ci 
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(solid line in Fig. 2(b) ) and the frequency line (dot line 
in Fig. 2(b) I. For frequencies far below the plasma fre- 



quency of the metal (e.g. below 100 THz), 62 can be ap- 



proximated as €2 = eo( 



With this approxima- 



tion, Lc becomes the well-known kinetic inductance L, 

2r92 19 ^ 



j^^^^a , which is non-dispersive as shown in Fig. 2(a) 

(dash line) . The resonance frequency can then be directly 
expressed as l/-\/(Lg + Lc)Ci (dash line in Fig. |2(b)| , 
and the quality factor Q — R^^y/{Lg + Lc)/Ci. For the 
2D SRR as shown in Fig. [TJ when we scale all the ge- 
ometrical parameters (a,6, r), Lc(equal to L^), Ci and 
R keep constant, while Lg scales proportional to the 
scale factor as a^. For silver SRR structure operating 
at optical frequencies, is comparable with or even 
larger than Lg, thus reducing Lg by size scaling is not 
effective to push resonant wavelength to short optical 
wavelength. For frequencies not far from the plasma 
frequency of the metal, Lc increases quickly when the 
wavelength decreases (cf. solid line in Fig. 2(a) I, which 
makes it even more difficult to push resonant wavelength 
to shorter optical wavelength by size scaling. We can see 
from Fig. 2(b) by comparing the solid line and dash line 



that for the given structure parameters, the resonance 
frequency determined by LC model with the approxima- 
tion Lc = Lc is over-estimated by about 50 THz. There 
have been efforts to reduce Lc by using metal with larger 
plasma frequency, e.g. aluminumii. However, the loss of 
aluminum is much larger than silver, which means that 
though magnetic resonance at short wavelength can be 
achieved, it is very weak. For silver SRR, another prob- 
lem during scaling is that although the quality factor 
decrease slowly, F' decrease very quickly, which reduces 
the amplitude of the magnetic resonance. So the con- 
ventional metallic SRRs have their limitations for appli- 
cations at short optical wavelength. In fact, there are 
two more factors that further reduce the amplitude of 
the magnetic resonance at short optical wavelength. In 
our analysis above, there are two assumptions that are 
less appropriate when |e2| is not large enough (e.g. when 
approaching the plasma frequency) as compared with the 
permittivity of the background material. First, the elec- 
tric field would not be strongly confined inside the ring. 
Second, apart from the displacement current inside the 
ring, the displacement current near the ring, whose di- 
rection is opposite to the displacement current inside the 
ring, should also be added to the total current in Eq. [TJ 

In Fig. [31 we show the scaling process of the SRR struc- 
tures in term of magnetic resonance frequency, where the 
solid lines represent simulation results while the dash and 
dot lines represent LC model results. From the simula- 
tion results, we can see the saturation phenomenon when 
the scaling factor a is small. SRRs with 9i ~ 20° sat- 
urates at about 500 THz, whereas SRRs with 9i = 40° 
saturates at about 600 THz. This is because 9i gives 
smaller Ci. We found again in this figure the different 
results given by LC model with (black dash line) and 
without (black dot line) the approximation Lc = Lc. 
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FIG. 3. (color online). The scaling of the simulated (solid 
lines) and LC model calculated (non-solid lines) magnetic res- 
onance frequency of conventional SRRs as a function of the 
unit cell size a. The black and red solid lines are simulated re- 
sults for SRRs with 9i = 20° and 40° , respectively. The dash 
line is the LC model calculated result for SRRs with 61 — 20° . 
The dot line is the LC model calculated result with the ap- 
proximation Lc = Lg for SRRs with 9i = 20°. The dash-dot 
line is the LC model calculated result with Ci doubled for 
SRRs with 01 = 20°. 



When the scaling factor 1/a is small (e.g. < 6/im~^), 
the over-estimation is not apparent, which means that 
in this region dispersive Lc can be well reduced to non- 
dispersive Lc. However, when the scaling factor is larger, 
the over-estimation is quite apparent, up to 130 THz. In 
this region, we have to use dispersive inductance. We 
would like to point out that the LC model (black dash 
line) can well describe the saturation phenomenon, al- 
though there is a frequency shift compared to the sim- 
ulation result (black solid line). Such a frequency shift 
between our LC model and simulation is not surprising. 
We have assumed a uniform distribution in the radial di- 
rection of the electric field inside the ring to get Eq. [1] 
whereas the numerical simulation shows that the electric 
field is far from uniform. This affects the estimation of 
the effective capacitance Ci in our model, as has been 
evidenced by Ref. [l^ The frequency shift is larger for 
larger 61 . Keeping this in mind, for SRRs with 9i = 20° 
we can simply double our effective capacitance Lc to well 
reproduce (black dash-dot line) the simulation result. In 
fact, we do not intend a very accurate theoretic model 
which requires lots of detailed factors to be involved at 
the cost of losing simplicity and physical intuitiveness to 
some degree. Our LC model aims at revealing the reso- 
nance conditions and key infiuential factors thus inspir- 
ing new designs with specific features of performance, e.g. 
magnetic metamaterial for operation at short wavelength 
as described in Sec. IIVI 
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FIG. 4. (Color online), (a) The scaling of the simulated 
(solid lines) and LC model calculated (dash lines) magnetic 
resonance frequency of HRRs as a function of the unit cell size 
a; (b) The effective permeability (retrieved from simulated S 
parameters) of HRRs with a — 90nm (denoted by the vertical 
line in panel (a)) and 6i = 40°, 20°, 10°, 5°, from left to right. 



The effective permeability has the same form as Eq. |4] 



fJ'eS = 1 



(5) 



where F' 



-F and R 



2r( 



i^j. For frequencies far 



below ujp, we can approximate ei as eo(— ^^j^*.^^ ) so that 

Lc — -^^^ and the resonant frequency can be expressed 

as l/-\/ (Lg + Lc)C2- We can now analyze why HRR is 
able to operate at short wavelength with strong magnetic 
response. First, is proportional to the length of the 
metal part, i.e. ^i, hence Lc of HRR is much smaller 
than that of the conventional SRR and we can thus ex- 
pect a much shorter resonant wavelength. Second, we can 
expect a higher quality factor of resonance because the 
resistance R is also proportional to 9i. Third, a smaller 
Lc also gives a larger F', thus a larger amplitude of mag- 
netic resonance. Additionally, unlike conventional SRR, 
HRR doesn't become dysfunctional when the permittiv- 
ity of the metal is too small to confine electric field near 
plasma frequency. This is because the high index di- 
electric material plays the dominant role in confining the 
electric filed inside the ring and the displacement current 
near the ring has the same direction as the displacement 
current inside the ring. 

The scaling of HRR is represented in Fig. 
6»i = 40° (black lines), 20° (red lines), 10° (greeFTTnes) 
and 5° (blue lines) by using LC model (dash lines) and 
simulation (solid lines) results. For each scaling curve 
(i.e. constant 0i), we still observe the saturation phe- 
nomenon as expected by our LC model. However, the 
saturation frequencies are much larger than those of con- 
ventional SRRs (cf. Fig. |3]) mainly owing to the reduc- 
tion of Lc. If we consider 6i as a variable parameter, we 
can see that the HRR is able to overcome the saturation 
problem. To demonstrate this ability, we keep a constant 
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IV. HYBRID METAL-DIELECTRIC RING 
RESONATOR TO OVERCOME THE 
SATURATION PROBLEM 

In order to operate at short wavelength, we proposed 
a new type of ring resonator, of which £2 is high index 
dielectric material while ei is metal. This is the case 
ei < 0, e2 > 0. We refer to this new type of ring resonator 
as HRR (Hybrid Ring Resonator) for short in this paper. 

For HRR, ei = eo(eoo — 1^2 j^^^^ )- For demonstration, 
we assume £2 = 6eo, which is reasonable since low- loss 
dielectric materials with permittivity around this value at 
short wavelength are available by using e.g. Ti02, SiC or 
ZnS from various deposition technologies. The effective 
capacitance C2 — (6e2)/(2r6'2) is positive while C2 ~ 
(6ei)/(2r6'i) is negative and is equivalent to a positive 
inductance Lc in series connection with a resistance R. 
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FIG. 5. (Color online). The simulated effective permeability 
of HRR using experimental material data of Ag from litera- 
ture. 
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as 90 nm (as indicated by the vertical line in Fig. 4(a) ) to higher loss for the experimental Ag above 800 THz. 



and decrease 9i from 40° to 5°. The effective permeabil- 
ity spectra are plotted in Fig. |4(b)[ By changing 9i , the 
resonance frequency increase from 800 THz to 1050 THz. 
It's worthwhile to point out that according to our Drude 
model of Ag, the real part of the permittivity is only -0.5 
at 1050 THz and is positive at larger than 1115 THz. We 
can almost say that we managed to push magnetic reso- 
nance frequency to as high as the permittivity of metal 
is negative. Additionally, as expected by our analysis 
with LC model, the amplitude and quality factor of the 
magnetic resonances keep quite healthy, only decreasing 
gently while being pushed to higher frequency. 

Till now we have simulated HRRs using Drude model 
in order to compare with the LC model and conventional 
SRRs. However, we know that due to the interband 
transitions, the fitted Drude model is no longer accurate 
above 700 THz. Therefore, we check again the perfor- 
mance of HRR by simulation using experimental mate- 
rial data of Ag2S,. There are two main differences between 
the experimental data and the Drude model. First, the 
frequency point where the permittivity change from neg- 
ative to positive is about 920 THz for the experimental 
data, which is much lower than Drude model (1115 THz); 
Second, the loss is several times higher at frequencies 
larger than 800 THz for the experimental data. The sim- 
ulated effective permeability spectra are shown in Fig. [S] 
The magnetic resonance can still be pushed to frequency 
close to 920 THz, although the amplitude and quality fac- 
tor are not good as the case using Drude model Ag due 



V. 



CONCLUSION 



We have proposed a model for ring resonators, based 
on which we summarized three cases of magnetic res- 
onances and studied in detail two cases in optical fre- 
quency. In our LC model, we utilized the effective kinetic 
inductance Lc, which is more accurate than L^. in describ- 
ing the performance of ring resonators, especially at short 
wavelength. In the framework of our model, we revisited 
the scaling of conventional SRRs and proposed a hybrid 
metal-dielectric ring resonator which is able to overcome 
the saturation problem of SRRs hence a very elegant de- 
sign to operate at short wavelength, down to ultraviolet 
wavelength. Our model is physically intuitive to guide 
magnetic metamaterial design in optical frequency. Our 
proposal of hybrid ring resonator design is just a proto- 
type, variations and improvements are possible for better 
performance and ease of fabrication. 
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